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Abstract 

After a pedagogical review of the simple constituent quark model 
and deep inelastic sum rules, we describe how a quark sea as produced 
by the emission of internal Goldstone bosons by the valence quarks can 
account for the observed features of proton spin and flavor structures. 
Some issues concerning the strange quark content of the nucleon are 
also discussed. 
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We shall first recall the contrasting concepts of current quarks vs con- 
stituent quarks. In the first Section we also briefly review the successes and 
inadequacies of the simple constituent quark model (sQM) which attempts 
to describe the properties of light hadrons as a composite systems of u, d, and 
s valence quarks. Some of the more prominent features, gleaned from the 
mass and spin systematics, are discussed. In the Sec. 2 we shall provide a 
pedagogical review of the deep inelastic sum rules that can be derived by way 
of operator product expansion and/or the simple parton model. We show in 
particular how some the sum rules in the second category can be interpreted 
as giving information of the nucleon quark sea. In the remainder of these 
lectures we shall show that the account of the quark sea as given by the chiral 
quark model is in broad agreement with the experimental observation. 



1 Strong Interaction Symmetries and the Quark 
Model 

In the approximation of neglecting the light quark masses, the QCD La- 
grangian has the global SU (3) L x SU (3) R symmetry. Namely, it is invariant 
under independent SU (3) transformation of the three left-handed and right- 
handed light quark fields. This symmetry is realized in the Nambu-Goldstone 
mode with the ground state being symmetric only with respect to the vector 
SU (3) L+iJ transformations. This gives rise to an octet of Goldstone bosons, 
which are identified with the low lying pseudoscalar mesons (it, K, rj) . For 
a pedagogical review see, for example, ref.[|l|. 

1.1 Current quark mass ratios as deduced from pseu- 
doscalar meson masses 

The light current quark masses are the chiral symmetry breaking parameters 
of the QCD Lagrangian. Their relative magnitude can be deduced from the 
soft meson theorems for the pseudoscalar meson masses. 

The matrix element of an axial vector current operator A 1 taken between 
the vacuum and one meson <f) b state (with momentum k^) defines the decay 
constant f a as 







4" 

A 4 



(j) (k)) = ik^fadat 
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where the SU(3) indices a, b.. .range from 1,2, ....8. This means that the di- 
vergence of the axial vector current has matrix element of 







<f) b (k)) = m 2 J a 5 ab , 



(1) 



If the axial divergences are good interpolating fields for the pseudoscalar 
mesons, we have the result of PCAC: 



&AI = m 2 J a 4> a . 



(2) 



Using PCAC and the reduction formula we can derive a soft-meson theorem 
for the pseudoscalar meson masses: 



mlfaSab 



-I 



- 



d xe 



-ik-x 







S(x ) ^(x),9^(0) 







Q 5b , Q 5a ,H(0) 



5a 







(3) 



where the axial charge is related to the time component of the axial vector 
current as Q 5a = J d 3 xA% (x) . 

If we neglect the electromagnetic radiative correction, only the quark 

masses, 

TC m = m u itu + mddd + m s ss, (4) 

break the chiral symmetry. Hence only such terms are relevant in the com- 
putation of the above commutators. [In actual computation it is simpler if 
one takes Ti m and Q 5a to be 3 x 3 matrices and compute directly the anti- 
commutator in <?4r7 7 5 (/, q\ b q = — \q |A a , A 6 | 7 5 g.] In this way we obtain: 

flml = ^ (m u + m d ) (o | (uu + dd) o) 

p2 2 



2 

- (m u + m s ) (0 \(uu + ss)\ 0) 



(5) 



2_2 
V 
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TO,, 



+ m d ) (o (uu + ddj 



-m s (0 \ss\ 0) 



1.1.1 Gell-Mann-Okubo mass relation and the strange to non- 
strange quark mass ratio 

Since the flavor SU(3) symmetry is not spontaneously broken, 



(0 \uu\ 0) = (0 dd 0) = (0 \ss\ 0) = /i 



(6) 



and f w = f K = L = /; Eq.QEJ) is simplified to 



m,. 



2m. 



"P 



(m n + m s ) 

2 u 3 
- (m n + 2m s ) — 



(7) 



m„. From this, we 



where we have made the approximation of m u ~ m^ 
can deduce the Gell-Mann-Okubo mass relation for the 0~ mesons: 



3m„ = 4m, 



m„ 



as well as the strange to nonstrange quark mass ratio 

m n _ m u + m d _ m 2 ^ J_ 
m, 2m, 2mi- — m£ 25 



(8) 



(9) 



1.1.2 Isospin breaking by the strong interaction & the m^/m^ ratio 

In order to study the ratio oim u /md, we need to include the electromagnetic 
radiative contribution to the masses. The effective Hamiltonian due to virtual 
photon exchange is given by 



h 7 = e 2 J d 4 xT ( j; m (x) J e v m (0)) 



(10) 



where D^ v (x) is the photon propagator. Thus, beside the contribution from 
Hm, we also have the additional term on the RHS of eqn.(|3|): 



a 



ab 







5a 







Now we make the observation (Dashen's theorem^) 
neutral mesons, we have [Q 5ct ,7i 7 ] = 0, which leads to 



(T 7 | TC 



(vr ) = a, (K°) = a, (rj) = 0. 
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For the electrically 
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On the other hand, J^ m is invariant (i.e. U-spin symmetric) under the in- 
terchange d <-> s, which transforms charged mesons tt + and K + into each 
other: 



(J 7 | 7T 



(5 " 



K~ 



(13) 



Consequently, we obtain the generalization of (|7]) as 
/ 2 m 2 (7T+) = (m u + m d ) /i 3 + /i 3 



f 2 m 2 (n°) = (m u + m d )fj, 3 
fm 2 (K + ) = (m u + m s )^ + fl 3 1 (14) 
f 2 m 2 (K°) = (m d + m s )fi 3 



f 2 m 2 {r\) = — {m u + nid + 4m s ) /i 3 . 

From this we can obtain the current quark mass ratios: 

m 2 (K°) + [m 2 (K+) - m 2 (tt+)] m s 
m 2 (K°) - [m 2 (K+) - m 2 (vr+)] ~ ma 

m 2 (K°) - [m 2 (K+) - m 2 (tt+)] txh 
[2m 2 (tt°) - m 2 (K )] + [m 2 (K+) - m 2 (tt+)] m u 

If we assume, for example, m s ~ 190 MeV, these ratios yield: 

m u ~5.3MeK m s ~9.5MeV or m n = 7 A MeV, (17) 

which are indeed very small on the intrinsic scale of QCD. This explains why 
the chiral SU(2) and isospin symmetries are such good approximations of the 
strong interaction. 



20.1 (15) 
1.8 (16) 



1.2 Quark masses from fitting baryon masses 

For the baryon mass we need to study the matrix elements {B \7H\B) . The 
flavor SU (3) symmetry breaking being given by the quark masses (|4]), we 
need to evaluate the matrix elements of the quark scalar densities u a between 
baryon states : 

lH. m = m u uu + m d dd + m s ss 
= m u + m 3 u s + m 8 ug 

uo = uu + dd + ss 

u 3 = uu — dd (18) 
u 8 = uu + dd — 2ss 



with 

m = | (m u + m d + m s ) 

m 3 = ~ (m u - m d ) 
m 8 = \ (m u + m d - 2m s ) 
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where, instead of the standard u a = q\ a q (A a being the familiar Gell-Mann 
matrices), we have used, for our purpose, the more convenient definitions 
of scalar densities by moving some numerical factors into the quark mass 
combinations mo,3,8- 

We shall first concentrate on the low lying baryon octet which, being the 
adjoint representation of SU (3) , can be written as a 3 x 3 matrix 



( ,/Je° + J\k s+ p \ 



B 



(19) 



The octet scalar densities u a can be related to two parameters (Wigner-Eckart 
theorem): 

(B \u a \ B) — atr (B ] ii a B) + (3tr (b ] Bu^j 

where u a is the scalar density expressed as a 3 x 3 matrix in the quark 
flavor space. The linear combinations (a ± (3) /2 are the familiar D and F 
coefficients. For example, we can easily compute: 

(p\u 8 \p) = a -2f3=(3F-D) mass (20) 
(p\u 3 \p) = a = (F + D) mass . (21) 

In this way the baryon masses with their electromagnetic self-energy sub- 
tracted (as denoted by the baryon names) can be expressed in terms of three 
parameters 

p = M + (a - 2/3)m 8 + am 3 (22) 

n = M.q + (a — 2/9)m 8 — am 3 
£+ = M + (a + (3)m 8 + {a - /?) m 3 

S° = M +(a + [3)m 8 

E~ = M + (a + /3)m 8 - (a - (5) m 3 

ET = Mo + (/3 - 2a)m 8 + f3m 3 

S° = Mo + (P - 2a)m 8 - (3m 3 

A = M -{a + (3)m 8 

We have 8 baryon masses and three unknown parameters Mo, ol and (5 - 
hence 5 relations, one of them should yield quark mass ratio m 8 /m 3 . 
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The ("improved") Gell-Mann-Okubo mass relation 



n + S~ = - (3A + 2S+ - E°) (23) 



2 

The Coleman-Glashow (U-spin) relation 

S" - H° = (p - n) + (E- - S+) (24) 

Absence of isospin 1 = 2 correction (i.e. M3 being a member of J = 1): 

5T + £ + - 2S° = (25) 

The hybrid relation: 

p — n 2~ — 5° . 

It should not be surprising that we have a relation relating S'C/ (2) 
breakings to SU (3) breakings, since M3 and us belong to the same 
octet representation. Recall that here the electromagnetic contribu- 
tion must be subtracted from our masses (sometimes called the tadpole 
masses). Since there is no Dashen theorem for the electromagnetic con- 
tributions to baryon masses, we must resort to detailed (& less reliable) 
model calculations. We quote one such result [Q for the electromagnetic 
contributions (AM) : 

p — n = (p — n) obs — (p — n) ~ —1.3 — 1.1 ~ — 2.4 Me V 



6.4- 1.3 ~ 5.1 MeV 



obs V / 7 

which yields ~ 0.02 on both sides of Eqn.(|2 

Both sides of Eq.(^) are related to the quark mass ratio 27713/ ~~ m 3) • Thus 
the above result leads to 

m u — rrid\ 

— ~ 0.02 27 

- m d - m s J B 

which is compatible with the current quark ratio deduced from pseu- 
doscalar meson masses Eqs.(|l5|) and fllB]) : 

' m u - m d \ ts - 1 
.m d - mj ps 



'-) ~ ~ 0.023. (28) 

/„, 1-20.1 v ; 
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1.3 The constituent quark model 

1.3.1 Spin-dependent contributions to baryon masses 

The sQM which attempts to describe the properties of light hadrons as a 
composite systems of u, d, and s valence quarks. The mass relations derived 
above may be interpreted simply as reflecting the hadrons masses as sum of 
the corresponding valence quark masses. For a general baryon, we have 

B = Mo + M x + M 2 + M 3 (29) 

where Aio is some SU (3) symmetric binding contribution. M\^^ are the con- 
stituent masses of the three valence quarks. We shall ignore isospin breaking 
effects: M u = Md = M n , and write the octet baryon masses as, 

N = M + 3M n (30) 

A = M + 2M n + M s 

E = M + 2M n + M s 

E = M + M n + 2M S , 

and the decuplet baryon masses as, 

A = M + 3M n (31) 

S* = M + 2M n + M s 

~* = M Q + M n + 2M S 

n = M + SM S . 

While it reproduces the GMO mass relations respectively, for the octet: 

iV + H = i(3A + S), (32) 

and for the decuplet (the equal-spacing rule): 

A - S* = S* - E* = E* - n, (33) 

it also leads to a phenomenologically incorrect result of A = S (reflecting 
their identical quark contents). Similarly, such a naive picture would lead us 
to expect that the N, E, E baryons having comparable masses as A, E*, E*. 
Observationally the spin 3/2 decuplet has significantly higher masses than 
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the spin 1/2 octet baryons. Similar pattern has also been observed in the 
meson spectrum: the spin 1 meson octet is seen to be significantly heavier 
than the spin mesons: M p> k*,u ^ M^^k^ even though they have the same 
quark contents. This suggests that there must be important spin-dependent 
contributions to these light hadron masses]^]. We then generalize Eq. (|29"D to 



B = Mo + Mi + M 2 + M 3 + k 



si-sa 



+ 



/ s 3 -s 2 



M X M 2 ) VM 3 M 



SrS3 

MiM 3 



(34) 



where is the spin of i-th quark, and the constant k one would adjust to 
fit the experimental data. This spin dependent contribution is modeled after 
the hyperfine splitting of atomic physics. For hydrogen atom we have a two 
body system hence only one pair of spin-spin interaction: M x = m e and 
M 2 = M p . The J arises from /x • B ~ /^ e -/x p /r 3 with the proportional 

constant worked out to be 



l^(0)| S 



where /j, = 2.79 is the magnetic moment of the proton in unit of nucleon mag- 
neton, and ip (0) is the hydrogen wave function at origin. Such an interaction 
accounts for the 1420 MHz splitting between the two IS states, which gives 
rise to the famous 21 cm line of hydrogen. For the case of baryon, one usually 
attributes such interaction to one-gluon exchange; but we shall comment on 
this point in later part of these lectures, at the end of Sec. 3.2. 
To compute the m*m- terms we need to distinguish three cases: 



The equal mass case: Mi = M 2 = M3 = M 



Si-S 2 



/ s r s 2 



M\Mi) VMiM 2 



+ 



1 



2M 2 



srgg 
MiM 2 

1 




2M 2 



[S(S+l)-3s(> + l) 



4Af 2 
3 



4Af 2 



for 
for 



S 
S 



1/2 
3/2 



(35) 



This is applicable for the N, A, Q baryons. 

(b) The unequal mass case, for example, (ssn) : Because of color an- 
tisymmetrization, the baryon wavefunction must be symmetric under the 
combined interchange of flavor and spin labels. Since we have a symmetric 
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superposition of flavor states, the subsystem (ss) must have spin 1, Namely, 
s s • s s = |(2 - 2s 2 s ) = J, and 



2s s ■ s n 



5>. -S 



-§-! = -! for 5=1/2 
= +| for S = 3/2 



or 



/ s l' s 2 



VMiM 2 
• 1 



s l' s 2 



+ 



/ Si-S2 



MiM 2 ) VMiM 2 
W ~ lo1 ' ' S = '/- 



+ 



l 



l 

4Af? 1 2M a M„ 



for 5 = 3/2 



(36) 



This case is applicable to H and H*, as well as E and E* because the sigma 
baryons are isospin 1 = 1 states (hence symmetric in the nonstrange flavor 
space) . 

(c) The A baryon: Because A is an isoscalar, the subsystem must be in 
spin state. From this one can easily work out the spin factor to be —4^2, 
independent of M s . 

Putting all this together into Eq . (|34]) we obtain, for the octet baryons: 



N 

A 
E 



M + 3M n 



3k 
~AMl 



M + 2M n + M s 
Mo + 2M n + M s + 



3 k 
AM? 



M + M n + 2M S + 



AMI M s Mn 



4M„ 2 M s Mn 



(37) 



and for the decuplet baryons: 



A 

E* 



Mo + 3M n + 



3k 
4MZ 



M + 2M n + M S + 
M + M n + 2M S + 



+ 



AMI 2M s M n 

K K 

AM] + 2M s M n 



n 



Mo + 3M S + 



3k 
AM?' 



(38) 
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One can obtain an excellent fit (within 1%) to all the masses with the pa- 
rameter values (e.g. 0) M.® — 0, jfe = 50 MeV and the constituent quark 
mass values of 

M n = 363 MeV, M s = 538 MeV. (39) 

Similarly good fit can also be obtained for mesons, with an enhanced 
value of k. Besides some different coupling factors this may reflect a larger 
|?/>(0)| 2 oc -R" 3 , which is compatible with the observed root mean square 
charge radii of mesons vs baryons: R m eson — 0.6 fm vs Rbaryon — 0.8 fm. 

1.3.2 Spin and magnetic moments of the baryon 

Another useful tool to study hadron structure is the magnetic moment of the 
baryon. Their deviation from the Dirac moments values (es/2M^) indicates 
the presence of structure. In the quark model the simplest possibility is that 
the baryon magnetic moment is simply the sum of its constituent quark's 
Dirac moments. Clearly, the magnetic moments are intimately connected to 
the spin structure of the hadron. Hence, we shall first make a detour into a 
discussion of the baryon spin structure in the constituent quark model. 



Quark contributions to the proton spin Because it is antisymmetric 
under the interchange of quark color indices, the baryon wavefunction must 
be symmetric in the spin-flavor space. Mathematically, we say that the 
baryon wavefunction should be invariant under the permutation group S3 - 
the group of permuting three quarks with spin and isospin labels. 

We shall concentrate on the case of proton. While the product wavefunc- 
tion is symmetric, the individual spin and isospin wavefunctions are of the 
mixed-symmetry type. There are two mixed-symmetry spin 
combinations: 



I wavefunction 



(i) x,s — symmetric in the first two quarks: Namely, the first two 
quarks form a spin 1 subsystem: (Notation for the spin-up and 



-down states: 



i, +7;) = a and |,— 1)=/3) 



|l,+l) = a 1 a 2 , |1,0) = —j= («i/3 2 +/3 1 a 2 ) 



|l,-l)=/3 1 /3 2 

which is combined with the 3rd quark to form a spin | proton: 



1 !\ 




1 1 

2' _ 2 




1 1 

2' + 2 
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or 

Xs = -^(2ai«2/?3 - ®iP 2 ®3 - Pia 2 a 3 ). (40) 

(ii) Xa " " antisymmetric in the first two quarks: The first two 
quarks form a spin subsystem: 



i + i).r |0 ' 0) 



1 1 

2' + 2 



or 

= ^( a iP 2 -/9ia 2 )« 3 - (41) 

While Xsa are the spin-\ wavef unctions, with identical steps, we can 
construct the two mixed-symmetry isospin-^ wavefunctions x's,A '■ 

Xs = -^(2MiM 2 4 - Uid 2 U 3 - rfiM 2 M 3 ) 

Xa = ~~^ Uld2 ~ rf l M 2)«3- (42) 

Both the spin wavefunctions (xs Xa) an d ^ ne isospin wavefunctions (x's Xa) 
form a two dimensional representation of the permutation group S3. For 
example, under the permutation operations of P\ 2 and P± 3 



a -( ; -. ) CD «•(::) -( 4 ".T ) (i 



M12 Ml 3 

where Mjj are 2-dimensional representations in terms of orthogonal matrices. 
Consequently, we find that the combinations such as (x% + x\) > (x's + Xa) 
and (XsXs + XaXa) are invariant under 5*3 transformations. In this way we 
find the symmetric proton spin-isospin wavefunction: 

\P+) = ^(XsX's + XaXa) (43) 
1 1 

= —^[-(2a 1 a 2 P 3 - a^ 2 a 3 - P 1 a 2 a 3 )(2u 1 u 2 d 3 - Wid 2 w 3 - d x u 2 u 3 ) 
1 

+ 2( a 1^2«3 - (3lCi 2 Cl?,){uid 2 U 3 - diU 2 U 3 )\ 
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[4 {u+u+d- + u + d^u + + d-U + u + ) 



QV2 

—2{u + U-d + + U-d + u + + d + u + u~ 
+u^u + d + + u + d + u^ + d + U-U + )} 

where we have used the notation of au = u + , (3d = gL, etc. In calculating 
physical quantities, many terms, e.g. tt+w+oL, u + d_u + and d_u + u + yield 
the same contribution. Hence we can use the simplified wavefunction: 

\p + ) = —j={2u + u + d_ — u + u_d + — u_u + d + ) (44) 
v6 

From this we can count the number of quark flavors with spin parallel or 
antiparallel to the proton spin: 

u+ = -, = -, d + = -, u_ = - (45) 

summing up to two u and one d quarks. From the difference 

Aq = q + - q„ (46) 

we also obtain the contribution by each of the quark flavors to the proton 
spin: 

4 1 

Au = - Ad= — As = 0, and AS = 1, (47) 
3 3 

where AS = Am + Ad + As is the sum of quark polarizations. 

Quark contributions to the baryon magnetic moments Instead of 
proceeding directly to the results of quark model calculation of the baryon 
magnetic moments, we shall first set up a more general framework. This will 
be useful when we consider the contribution from the quark sea in the later 
part of these lectures. We shall pay special attention to the contribution by 
antiquarks. If there are antiquarks in the proton, the definition in Eq.(PB]) 
becomes 

A<7= (q+ + (q+-q-) =A q + A w (48) 

Thus the quark spin contribution Aq is the sum of the quark and antiquark 
polarizations. For the g-flavor quark contribution to the proton magnetic 
moment, we have however 

(q) = A q H q + Ag/Mg = (A q - A ? ) fi q = Aqfl q (49) 



14 



where fi q is the magnetic moment of the g-flavor quark. The negative sign 
simply reflects the opposite quark and antiquark moments, /% = —fi q . Thus 
the spin factor that enters into the expression for the magnetic moment is 
Aq, the difference of the quark and antiquark polarizations. If we assume 
that the proton magnetic moment is entirely built up from the light quarks 
inside it, we have 

/i p = Aufi u + Adfi d + As/i s . (50) 

In such an expression there is a separation of the intrinsic quark magnetic 
moments and the spin wavef unctions. Flavor-S'[/(3) symmetry then implies, 
the proton wavefunction being related the S + wavefunction by the inter- 
change of d <-> s and d ^ s quarks, the relations (^Au^j^ + = (Au^J = Au, 

(^Adj^ + = As, and (As)^ + = Ad; similarly it being related to the S° wave- 
function by a further interchange of u <-> s quarks, thus ( Ac/) _ o = ( Ad^j 
As, (As)„ = (Au) s+ = Au, and (Au)„ o = (As) s+ = Ad. We have, 



E+ 



fM s + = Aufi u + Asfi d + Adfi s , (51) 
// s o = Adfi u + Asfi d + Aw/i s , (52) 



the intrinsic moments /i„ being unchanged when we go from Eq.(|50|) to 
Eqs. (|5l|) and (0). The n, E~, and S~ moments can be obtained from 
their isospin conjugate partners p, S + , and S° by the interchange of their 
respective u «-> d quarks: (^Auj^_ = (^Adj^ + = As, etc. 

fi n = Adfi u + Aufi d + Asfi s , (53) 
/x E _ = Asfi u + Aufi d + Ad[i s , (54) 
A*h- = Asfi u + Adfi d + Aufi s . (55) 

The relations for the I z = 0, K = moments are more complicated in 
appearance but the underlying arguments are the same. 

0a = ~(A^ + 4Ad + As)(ti u + Vd) (56) 
+- UAu - 2Ad + 4As) u, 

Hju = 0=(Au-2Ad + As)( f i u - f i d ). (57) 
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In the nonrelativistic constituent quark model, there is no quark sea and 
hence no antiquark polarization, Aq = 0. This means that in the sQM we 
have Aq = Aq. After plugging in the result of Eq. 
in the 2nd column of the Table 1 below: 



Baryon 


mag 
(* 


moment 
= » q ) 


it = -2d 
s = 2d/3 


d = —0.9fi N 


exptl # 


P 


(4m 


-d) /3 


-3d 


2.7 


2.79 


n 


(4d 


-u) /3 


2d 


-1.8 


-1.91 


£+ 


(4m 


-s)/3 


-26d/9 


2.6 


2.48 


£~ 


(Ad 


-s)/3 


10d/9 


-1.0 


-1.16 


CO 


(4s 


-u) /3 


14d/9 


-1.4 


-1.25 




(4s 


-d) /3 


5d/9 


-0.5 


-0.68 


A 




s 


2d/3 


-0.6 


-0.61 


AE 


(d- 


w)/V3 


\/3d 


-1.6 


-1.60 



Table 1 . Quark contribution to the octet baryon magnetic moments. 

Instead of trying to get the best fit at this stage, we shall simplify the 
result further with the following observation: Because of the assumption 
M u = Md, we have fi u = —2{i d . The proton and neutron moments are then 
reduced to « p = — 3u d and /i n = 2u d , and thus the ratio 

^ = -1.5 (58) 

which is very close to the experimental value of —1.48. Furthermore, we have 
seen in previous discussion that constituent strange-quark mass is about a 
third heavier than the nonstrange quarks M s /M n ~ 3/2, we can make the 
approximation of fi s = 2fi d /3. In this way, all the moments are expressed in 
terms of the d quark moment, as displayed in the 3rd column above. One 



47), we obtain the result 
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can then make a best over-all-fit to the experimental values by adjusting this 
last parameter fi d . The final results, in column 4, are obtained by taking 
fi d = —0.9/%, where /% is nucleon magneton e/2Mjv. They are compared, 
quite favorably, with the experimental values in the last column. We also 
note that, with the d quark having a third of the electronic charge, the fit- 
parameter of fi d = —0.9 fi N translates into a d quark constituent mass of 

= Mjv = 348 MeV and M = ^MlL = 522 MeV, (59) 
3 x 0.9 2 v ' 

which are entirely compatible with the constituent quark mass values in 
Eq.(^), obtained in fitting the baryon masses by including the spin-dependent 
contributions. 

1.3.3 sQM lacks a quark sea 

So far we have discussed the successes of the simple quark model. There 
are several instances which indicate that this model is too simple: sQM does 
not yield the correct nucleon matrix elements of the axial vector and scalar 
density operators. 

Axial vector current matrix elements The quark spin contribution to 
proton Aq in Eq.fl4*5|) is just the proton matrix element of the quark axial 
vector current operator 



2s M Ag = {p, s qi^l^q p, sj = 2s ^ (q + - g_ + q + - q_) (60) 

where s M is the spin-vector of the nucleon, as the axial current vector corre- 
sponds to the non-relativistic spin operator: 



<m 5 <? = ?M n - ( 61 ) 




Through SU (3) these matrix elements can be related to the axial vector 
coupling as measured in the octet baryon beta decays. In particular, we 
have 

(An-Ad) exptl = 1.26 
(Aw + Ad — 2As) expt = 0.6 (62) 
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which is to be compared to the sQM results of Eg. (fT7|) : 

(Au-Ad) sQM = 5/3 
Ad-2Asl.™, = 1. 



(Am 



>sQM 



(63) 



Scalar density matrix elements The matrix elements of scalar density 
operator qq can be interpreted as number counts of a guark flavor in proton 



(p \qq\p) = q + q 



(64) 



where q (q) on the RHS denotes the number of a guark (antiguark) flavor 
in a proton. Namely, the proton matrix element of the scalar operator qq 
measures the sum of guark and antiguark number in the proton (opposed to 
the difference q — q as measured by q^q ). It is useful to define the fraction 
of a quark-flavor in a proton as 



F(q) 



(p\qq\p) 


(p 


uu + dd + ss 


p) 



(65) 



We already have calculated proton matrix element of the scalar density in 
the subsection on the baryon masses, Egs.(|lD and ©. Thus we have 



^(3) 

F (8) F(u) + F (d) - 2F (s) 



F{u)-F (d) 



a: 



a: 



2(3 



(66) 



The parameters (a, (3) can be deduced from Eg.(|2^) in the SU (3) symmetric 
limit (msUs = 0) , as 



3m« 



P 



N 



3m» 



Thus the ratio 



£(3) 
F(8) 



M E -M s 



- exptl 



2M N -M s -M s 



0.23 



which is to be compared to the sQM value of 

F(3) 



F(8) 



J sQM 



(67) 



(68) 



The simplest interpretation of these failures is that the sQM lacks a guark sea. 
Hence the number counts of the guark flavors does not come out correctly. 
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1.4 The OZI rule 



The simple quark model of hadron structure discussed above ignores the 
presence of quark sea. Even when the issue of the quark sea in nonstrange 
hadrons is discussed, its (ss) component is usually assumed to be highly 
suppressed. This is based on the OZI-rule0, which was first deduced from 
meson mass spectra. In this Subsection we briefly review this topic. 

1.4.1 The OZI rule for mesons 

The three (qq) combinations that are diagonal in light-quark flavors are the 
two isospin 1 = 1 and states of a flavor-SU(3) octet together with a SU(3) 
singlet. Isospin being a good flavor symmetry, there should be very little 
mixing between the 1 = 1 and states. On the other hand, the flavor-SU(3) 
being not as a good symmetry, we anticipate some mixing between the octet- 
and the singlet- 1 = states: 



Pseudoscalar meson masses and mixings The Gell-Mann-Okubo mass 
relation for the 0~ mesons, before the identification of rj as the 8th member 
of the octet, may be interpreted as giving the mass of this 8th meson: 



which is much closer to the rj meson mass of m v = 547 Me V than = 
958 MeV. The small difference m 8 — can be attributed to a slight mixing 
between the octet and singlet isoscalars. Namely, we interpret rj and rf 
mesons as two orthogonal combinations of |8) and |0) with a mixing angle 
that can be determined as follows: 

/ m| mo 8 \ ( cos 9 sin 9 \ / \ / cos 9 — sin 9 \ 

I mg mQ J I — sin 9 cos 9 J I mi, J \ sin 9 cos 9 J ' 





(69) 




(70) 



Hence 




i.e. a small 9p ~ 11°. 



(71) 
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Vector meson masses and mixings We now apply the same calculation 
to the case of vector mesons: 



*2 

8 3 



\ (Ami, - m2 p ) = (929 MeVf 



which is to be compared to the observed isoscalar vector mesons of u (782 MeV) 
and (1020 MeV) . This implies a much more substantial mixing. The diag- 
onalization of the corresponding mass matrix: 



'80 

requires a mixing angle of 



mf mil ] , ( m l 

mil mf 



ml 



sin 2 e v = m l m \ or fl y ~50°. (72) 

mj - mi 

The physical states should then be 

\u) = cos0 y |8)+sin0y |0) |0) = - sin ^ y |8) + cos ^ y |0) . (73) 

After substituting in Eqs.(^) and ([72]) into Eq. (|73|) , we have 

|w) = 0.7045 \uu + dd) + 0.0857 \ss) \4>) = -0.06 \uu + dd) + 0.996 |ss) . 

This shows that w has little s quarks, while the <fi mesons is vector meson 
composed almost purely of s quarks. Such a combination is close to the 
situation of "ideal mixing" , corresponding to an angle of 6* — 55°, with the 
non-strange and strange quarks being completely separated: 

\u) = —j= uu + dd) \4>) = \ss) . (74) 
v2 ' 

The OZI rule It is observed experimentally that the 4> meson decay pre- 
dominantly into strange-quark-bearing final states, even though the phase 
space, with m^ > m u , favors its decay into nonstrange pions final states: 

uj ^3tt 89% -> KK 83% 

-> pTT 13% 

-> 3tt 3% 
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with a ratio of partial decay widths T(0 — > 37r)/T(c<j — > 37r) = 0.014. 

This property of the hadron decays has been suggested to imply a strong 
interaction regularity: the OZI-rule — the annihilations of the ss pair via 
strong interaction are suppressed J7j. We remark that this suppression should 
be interpreted as a suppression of the coupling strength rather than a phase 
space suppression due to the larger strange quark mass (i.e. it is above and 
beyond the conventional flavor SU(3) breaking effect.) 

The extension of the OZI-rule to heavy quarks of charm and bottom has 
been highly successful. For example it explains the extreme narrowness of the 
observed J/ ip width because this (cc) bound state is forbidden to decay into 
the OZI-allowed channel of DD because, with a mass of mm ~ 3100 Me V, 
it lies below the threshold of 2m£> ~ 3700 MeV. 

From the viewpoint of QCD, applications of the OZI-rule to the heavy 
c, b, and t quarks are much less controversial than those for strange quarks - 
even thought the rule was originally "discovered" in the processes involving 
s quarks. For heavy quarks, this can be understood in terms of perturbative 
QCD and asymptotic freedom^. It is not the case for the s quark which, 
as evidenced by the success of flavor-SU(3) symmetry, should be considered 
a light quark. Furthermore, the phenomenological applications of the OZI 
to strange quark processes have not been uniformly successful. In contrast 
to the case of vector mesons Eq.([72]), there is no corresponding success for 
the pseudoscalar mesons — as evidenced by the strong deviation from ideal 
mixing in the r\ and rf meson system, Eq. ([7T|) . 

1.4.2 The OZI rule and the strange quark content of the nucleon 

A straightforward application of the s quark OZI rule to the baryon is the 
statement that operators that are bilinear in strange quark fields should 
have a strongly suppressed matrix elements when taken between nonstrange 
hadron states such as the nucleon. In particular we expect the fraction of s 
quarks in a nucleon, Eq.(^), should be vanishingly small. 



s + s (N\ss\N) 



(75) 



E{q + q) (N\uu + dd + ss\N) 



The "measured" value of the pion-nucleon sigma term|9|]: 




(76) 
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and the SU(3) relation 



m s (N\u 8 \N} = -(m n - 
M k - M= ~ -200 MeV, 



N 



uu + dd — 2ss 



N 



(77) 



which is obtained from Eqs.(j20Q and fl22|) in the isospin invariant limit (1713U3 = 
allow us to make a phenomenological estimate of the strange quark content 
of the nucleon|l0|: We can rewrite the expression in Eq.(|75|) 



as 



F(s) 



(N 


(uu + ddj — 


uu + dd — 


2ss) 


N) 


(N 


3 [uu + dd^j — 


(uu + dd - 


- 2ss) 


n) 



a wN - 25 MeV 
3a nN - 25 MeV 



(78) 



where we have used (ITT]) and the current quark mass ratio mg/m s = —8 corre- 
sponding to m s /m n = 25 of Eq.(|9|). Thus the validity of OZI rule, F (s) = 0, 
would predict, through (1781), that a^N should have a value close to 25 MeV. 



However, the commonly accepted phenomenological value |T]J is more like 
45 MeV, which translates into a significant strange quark content in the nu- 
cleon: 

F(s)~0.18. (79) 

We should however keep in mind that this number is deduced by using flavor 
SU (3) symmetry. Hence the kinematical suppression effect of M s > M n has 
not been taken into account. 



2 Deep Inelastic Scatterings 

2.1 Polarized lepton-nucleon scatterings 

There is a large body of work on the topic of probing the proton spin structure 
through polarized deep inelastic scattering (DIS) of leptons on nucleon target. 
The reader can learn more details by starting from two excellent reviews of 



12 and 13 
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2.1.1 Kinematics and Bjorken scaling 



For a lepton (electron or muon) scattering off a nucleon target to produce 
some hadronic final state X, via the exchange of a photon (4-momentum 
q^), the inclusive cross section can be written as a product 



da (I + N -> I + X) oc l^Wf 



(80) 



where is the known leptonic part while W„v is the hadronic scattering 
amplitude squared, X)x \T (7* (?) + iV (p) —> X)\ , which is given, accord- 
ing to the optical theorem, by the imaginary part of the forward Compton 
amplitude: 



W, v = ±-lmf(p,s \T (j; m (x) J e u m (0)) I p, s) e^d'x 



-9iw + 



q 2 



Fi(q 2 ,v) 



p-q 



Pv 



p-q 

Q , 



F2 (q 2 , u) 



p-q 



B9i(q,v) g fj g 2 (q 2 ,v) 
+ p ■ qs p - s ■ qjf- -2- 

p-q (p-q) 



where 



-Q 2 < and v = 



p-q 

M 



(81) 



(82) 



M being the nucleon mass. s a = ujy (p, s) 7 a 7 5 Mjv ip, s) is the spin- vector of 
the proton, and the variable v is the energy loss of the lepton, v = E — E'. 
We have defined the spin-independent Ei^iq 2 ,^) and the spin-dependent 
9i,2(q 2 ,v) structure functions. In particular, the cross section asymmetry 
with the target nucleon spin being anti-parallel and parallel to the beam of 
longitudinally polarized leptons is given by the structure function gi : 



da n da" 



^ ME ln 

— xy (2-y) 9l + 



M 2 ^ 
M 2 , 



(83) 



where x = and y 
spin- asymmetry, 



dxdy dxdy nQ 4 
QL. cnri n, — ii . \ n practice one measures g\ via the (longitudinal) 

_ da^ 



2x^- 
2X F 2 



(84) 



do-W + da^ 

in the kinematic regime of v ^> y/Q 2 . 

To probe the nucleon structure at small distance scale we need to go 
to the large energy and momentum-transfer deep inelastic region — large 
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Q 2 and v, with fixed x. In the configuration space, this corresponds to the 
lightcone regime. The statement of Bjorken scaling is that, in this kinematic 
limit, the structure functions approach non-trivial functions of one variable: 

F lt2 (q 2 , v) -> F 1>2 (x) , g lt2 (q 2 , u) -> g lj2 (x) . (85) 

Such problems can be studied with the formal approach of operator product 
expansion, which has a firm field theoretical-foundation in QCD, or the more 
intuitive approach of parton model, which can lead to considerable insight 
about the hadronic structure. 

2.1.2 Inclusive sum rules via operator product expansion 

The forward Compton amplitude is the matrix element, taken between 
the nucleon states 

Tpv = {p,s\t^ u \p,s) , (86) 
of the time-order product of two electromagnetic current operators 

t iu , = ij d 4 xe^ x T ( J„ (x) J v (0)) . (87) 

It is useful to express the product of two operators at short distances as 
an infinite series of local operators, O a (x) Ob (0) = J2i Ci (x) Oi (0), as it 
is considerably simpler to work with the matrix elements of local operators 
Oi (0). For DIS study we are interested in the lightcone limit x 2 — > 0. Hence 
operators of all possible dimensions (dj) and spins (n) are to be included: 

O a (x) O b (0) = J2Q (x 2 ) O' 1 ->■■ (0) (88) 

i,n 

where (0) is understood to be a symmetric traceless tensor operator 

(corresponding to a spin n object). From dimension analysis we see that the 
coefficient 

q (x 2 ) ~ (y^) Ti - dA - dB 

where r« = d{ — n is the twist of the local operator ^ (0) . Thus in the 
lightcone limit x 2 — > 0, the most important contributions come from those 
operators with the lowest twist values. 

In the short distance scale, the QCD running coupling is small so that per- 
turbation theory is applicable. In this way the c-numbers coefficients Cj (x 2 ) 
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can be calculated with the local operators Of 1 " Mn (0) being the composite 
operators of the quark and gluon fields. 

We are interested, as in Eq.(^7]), in the operator products in the momen- 
tum space. Namely, the above discussion has to be Fourier transformed from 
configuration space into the momentum space: x — > q, with the relevant limit 
being Q 2 —>■ oo. The spin-dependent case corresponds to an operator product 
antisymmetric in the Lorentz indices \i and v : 



E C(3)(g 2 ,« s ) — 5 (89) 



i/>,n=l,3,. 



-Q 



where (7(3) (g 2 , a s ) = 1+0 (a s ) , [the subscript (3) reminds us of others terms, 
1 & 2, that contribute to the spin-independent amplitudes F\p\. O^""' 1 " is 
a twist- two pseudotensor operator: 

o%r" n = 4 (£T V ^ .... ^ n oo) 

where ^ is the quark field with charge e^. The crossing symmetry property 

t^{p,q) =t vli (p,-q) (91) 

implies that only odd-n terms appear in the \pu] series. (By the same to- 
ken, only even-n terms contribute to the spin-independent structure function 

Fl,2-) 

The spin-dependent part of the forward Compton amplitude Eq.flgSI) is 

T^y] = (p, s p, = ie^a^s ^ 1 q ,U + ... (92) 

Namely, lmgi(q 2 ,u) = 2ngi (q 2 , v) . When we sandwich the OPE terms 
Eqs.(|89"|) and ( PD| ) into the nucleon states we need to evaluate matrix ele- 
ment 



p,s 



n /3fM 2 ....^ n 

Art 



p,s) = 2e%A n ^p^....p^ (93) 



Plug Eqs.(P|) and <M) into Eq.(B3) we have 



ituva^s 13 -^- = V) ie ta/a pq a s p (p • q) n 1 2elA n ^ 

P-Q n=ll... \-9 J 
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or 

ft = X;2C(3)44^u; n (94) 

ifi,n 

where uj — _^ 2 is the inverse of the Bjorken-x variable. Asymptotic freedom 
of QCD has allowed us to express the structure function as a power series 
in o;,Eq. (|94) with calculable c number coefficients CV3) and "unknown" long 
distance quantities A n ^. To turn this into a useful relation we need to invert 
the summation over n {i.e. to isolate the coefficient A n ^). For this we can 
use the Cauchy's theorem for contour integration: 



i/^|lr = E2Q3)e^ ) (95) 

which can be related to physical processes by evaluating the LHS integral 
with a deformed contour so that it wraps around the two physical cuts, 
u = (1, 00) and (—00, —1) . (The second region corresponding to the cross- 
channel process.) Using 

gi (uj + ie) — g\{uj + ie) = 2i Im g\ {uj) = 4:ing 1 (uj) (96) 

and the crossing symmetry property 

9i (P, <?) = -9i (p, ~q) or g 1 (uj, g 2 ) = -g x (-uj, q 2 ) , (97) 

we then obtain 



2ni J u n+1 n Ji uj u+1 n J _oo uj u+1 

= 2 [1 - (-1)"] [°° duj^- = 4 f 1 x n ~ l 9l (x) dx. (98) 

Jl UJ n+1 Jo 

We recall that the spin-index n must be odd. The first-moment (n = 1) sum 

f dx 9l (x, Q 2 )=\Y. C ®A> A hi> (99) 
Jo 1 ^ 

is of particular interest because the corresponding matrix element on the 
RHS can be measured independently, Cf. Eqs.(|60|) and (|93|): 



2A l s fS = (p, s I^TsV'l P, s ) = (100) 
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Without including the higher order QCD corrections in the coefficient, we 
have the g\ sum rule for the electron proton scattering: 



£ dxg{ (x, Q 2 ) = \ (~Am + I Ad + iAs) . 



(101) 



For the difference between scatterings on the proton and the neutron targets, 
we can use the isospin relations (Au) = Ad and (Ad) = Am to get: 



dx 



g p l (x 1 Q 2 )-g n l (x,Q 2 )]=-(Au-Ad) 



(102) 



The matrix element on the RHS: 
2/ (Aw -Ad) = 



p, s M7 /3 7 5 m — d7 /3 7 5 d p, s 
(p, s M7 /3 7 5 d n, sj> = 2s p g A 



(103) 



is simply the axial vector decay constant of neutron beta decay. Including 
the higher order QCD correction to the OPE Wilson coefficient, one can then 
write down the Bjorken sum rule: 



dx 



(x,Q 2 )-< (x,Q 2 )] = 9 -^C iNS) 



with the non-singlet coefficient [14 



(TVS) 



l _ «s _ 43 /as 

7T 12 V 7T 



- 20.22 



7T 



+ ... 



(104) 



(105) 



All experimental data are consistent with this theoretical prediction. 

Remark Anomalous dimension and the Q 2 -dependence: The Q 2 -dependence 
of the moment integral, such as LHS of Eq. (|99|) , are given by a s (Q 2 ) ~ 
1/ In Q 2 in the coefficient function and by the Q 2 -evolution of the oper- 
ator according to the renormalization group equation [15], which yields 

(106) 



(p,s 


o\ Q 


P,s) 


' a s (Q) " 


2 

2b 


(p,s 


<?lQo 




a s (Q )_ 





where 7 is the anomalous dimension of the operator O and h is the 
leading coefficient in the QCD f3 function. The label Q in the matrix 



27 



elements refers to the mass scale at which the operator is renormalized, 
chosen at /i 2 ~ Q 2 in order to avoid large logarithms. For the gi sum 
rule Eq.(p9|) the Q 2 -dependence is particularly simple. The non-singlet 
axial current is (partially) conserved, hence has anomalous dimension 
7 = 0. The singlet current is not conserved because of axial anomaly 
(see discussion below). But it has very weak Q 2 -dependence because 
the corresponding anomalous dimension starts at the two-loop level. 



2.1.3 The parton model approach 



The gi sum rule of Eq.( |101| ) has been derived directly through OPE from 



QCD. We can also get this result by using the parton model, which pictures 
the target hadron, in the infinite momentum frame, as superposition of quark 
and gluon partons each carrying a fraction (x) of the hadron momentum. For 
the short distance processes one can calculate the reaction cross section as an 
incoherent sum over the rates for the elementary processes. Thus in Compton 
scattering, a photon (momentum g M ) strikes a parton (xp^) turning it into a 
final state parton (q^ + xp^), the initial and final partons must be on shell: 

2 2 ~Q 2 

(xPid =(qii + xPfi) or x = - . (107) 

Zp ■ q 

Hence the Bjorken-x variable has the interpretation as the fraction of the lon- 
gitudinal momentum carried by the parton. A simple calculation fl6|] shows 
the scaling structure functions being directly related to the density of partons 
with momentum fraction x : 

F£ (x) = x £ e 2 q [q(x) + q(x)} (108) 

q=u,d,s 

and 

gl ( x ) = \ 4 [?+ ( x ) - 1- ( x ) + <?+ (x) - q- {x)\ 

q=u,d,s 

= ^[^W + M^l^SW 1 ) (109) 

Thus the spin asymmetry of Eq.fl8~4|) has the interpretation as 

E q el[q(x) + q(x)] ■ (110) 
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Comparing this interpretation of the spin-dependent structure function 
to that for the proton matrix elements of the axial vector current Eq. (|60"D , 
we see that the g\ sum rule Eq. ( |101| ) implies the consistency condition of 

/ Q± ( x ) dx = q± / q± (x) dx = q±. (HI) 
Jo Jo 

In other words, the proton matrix element of the local axial vector current 
(p, s \OA,q \ P, s) can be evaluated, in the partonic language, by taking the axial 
vector current between quark states ((q, h \0 A , q \ q, h) = 2h) and multiplying 
it by the probability of finding the quark in the target proton: 

(p,s\0 A Jp,s) = Y / (<l,h\0 Atg \q,h)q h (x) = (Aq) p (112) 

q,h 

where (Aq) = Aq 

Aq (x) = q + (x) - g_ (x) + q + (x) - g_ (x) = A q (x) + A 5 (x) . (113) 



Ellis-Jaffe sum rule and the phenomenological values of Aq Besides 

Au - Ad = g A = F + D = 1.2573 ± 0.0028, (114) 

if we assume flavor SU(3) symmetry, we can fix another octet combination 

Au + Ad - 2As = A 8 = 3F - D = 0.601 ± 0.038 (115) 

which can be gotten by fitting the axial vector couplings of the hyperon beta 
decays]!?]]. In this way Eq. (|101j ) can be written as 



dxg\ [x) 



°^ (3 9A + A«) + %AS 



3(5 



9 



;ii6) 



where AS = Au + Ad + As. The non-singlet coefficient has been displayed 
in Eq.(|105|) while the singlet term has been calculated to be]|Hf| 



a 



(S) 



1 1.0959 — +... 



7T 



at* 
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(117) 



If one assume As = 0, thus AE = A 8 we then obtain the Ellis-Jaffe sum 
ra/e]p^] with the RHS of Eq.( |116j ) expected (for a s ~ 0.25) to be around 
0.175, had become the baseline of expectation for the spin-dependent DIS. 
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The announcement by EMC collaboration in the late 1980 's that it had 



extended the old SLAC result p0| to new kinematic region and obtained an 



experimental value for T p deviated significantly from the Ellis- Jaffe value pi] 
had stimulated a great deal of activity in this area of research. In particular 
another generation of polarized DIS on proton and neutron targets have been 



performed by SMC at CERN|22j and by E142-3 at SLAC0. The new data 
supported the original EMC findings of As 7^ and a much-less-than-unity 
of the total spin contribution AS <C 1, although the magnitude was not as 
small as first thought. The present experimental result may be summarized 
as [[24]] 



Am = 0.82 ±0.06, Ad = -0.44 ± 0.06, (118) 
As = -0.11 ±0.06, AS = 0.27 ±0.11. 



The deviation from the simple quark model prediction Eq.( 

(A?) exptl < (Ag) sQM (119) 

indicates a quark sea strongly polarized in the opposite direction from the 
proton spin. That the total quark contribution is small means that the proton 
spin is built up from other components such as orbital motion of the quarks 
and, if in the relevant region, gluons. 

2.1.4 Axial vector current and the axial anomaly 

The most widely discussed interpretation of the proton spin problem is the 
suggestion that the gluon may provide significant contribution via the axial 
anomaly P5|. Let us first review some elementary aspects of anomaly. The 



SU (3) coior gauge symmetry of QCD is of course anomaly-free. The anomaly 
under discussion is the one associated with the global axial U (1) symme- 
try. Namely, the SU (3)-singlet axial current = J2 q = u ,d,s Ql^l5Q nas an 
anomalous divergence 

& A f= £ 2m q (qi l5 q) + nArG^G fll/ (120) 

q=u,a,s 

where G^ u is the gluon field tensor, G^ u its dual, rif = 3 is the number of 
excited flavors. For our purpose it is more convenient to express this in terms 
of each flavor separately. 

0" (9Vy6?) = 2m, {qi lh q) + ^trCTG^ (121) 
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Axial anomaly enters into the discussion of partonic contributions to the 
proton spin as follows: Because anomaly, being related to the UV regular- 
ization of the triangle diagram, is a short-distance phenomena, it makes a 
hard, thus perturbatively calculable (though not the amount), contribution 
from the gluon so that Eq. (|112[) is modified: 



(p, s \0 A , q \p, s)=Y, (q, h \0 A>g \ q, h) Q h {x) +J2(G,h \0 A , q \ G, h) G h (x) 

q,h G,h 

(122) 

where Gh, just as the quark density Qh being given by Eq. (|113|) , is the spin- 
dependent gluonic density. The gluonic matrix element of the axial vector 
current (G,h\0 Atq \G,h) is just the anomaly triangle diagram which, with 
(q, h \ Atq \ q, h) normalized to ±1, yields a coefficient of ^f 4 . In this way the 
proton matrix element of the axial vector current is interpreted as being a 
sum of "true" quark spin contribution AQ and the gluon spin contribution: 

Aq (x) = AQ (x) - -^AG (x) , (123) 

where AG (x) = G + (x) — G_ (x) . Superficially, the second term is of higher 
order. But because the InQ 2 growth of AG (due to gluon blemsstrahlung 
by quarks) compensates for the running coupling a s ~ (InQ 2 ) 1 , the com- 
bination a s AG is independent of Q 2 at the leading order, and the gluonic 
contribution to the proton spin may not be negligible. However in order to 
obtain the simple quark model result of AS = 0, a very large AG is required: 

- —AG = As ~ -0.1 AG ~ 2.5. (124) 
2tt 

2.1.5 Semi-inclusive polarized DIS 

From the inclusive lepton nucleon scattering we are able to extract the quark 
contribution to the proton spin, Aq = A q + A q . Namely, we can only get the 
sum of the quark and antiquark contributions together. More detailed infor- 
mation of the spin structure can be obtained from polarized semi-inclusive 
DIS, where in addition to the scattered lepton some specific hadron h is also 
detected. 

l + N ^l + h + X 
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The (longitudinal) spin asymmetry of the inclusive process can be expressed 
in terms of quark distributions as in Eq. ( |110| ): 

A ^i^) (125) 

E g e 2 q {q + q) 

Similarly one can measure the spin-asymmetry measured in semi-inclusive 

case: 



where Dl 1 , the fragmentation function for a quark g to produce the hadron 
h, is assumed to be spin-independent. Separating A q from A q is possible 
because ^ D^. For example, given the quark contents such as n + ~ 
and tc~ ~ (w<i), we expect 

Df^Df, Df^Df, and « Df , Df < 



In this way the SMC collaboration |26| made a fit of their semi-inclusive data, 



in the approximation of A s = Aj and A s = A s oc s (x) (the strange quark 
distribution did not play an important role, and the final result is insensitive 
to variation of As). SMC was able to conclude that the polarization of the 
non-strange antiquarks is compatible with zero over the full range of x : 

A fi = A d - = -0.02 ± 0.09 ± 0.03 (127) 

This is to be compared to their result for Ag = A q — Ag : 

Au = 1.01 ± 0.19 ± 0.14 Ad = -0.57 ± 0.22 ± 0.11. 

Namely, while the data from inclusive processes suggest that the quark sea 
is strongly polarized — as indicated by the large deviation of measured Ag 
from their simple quark model prediction Eqs.( TTH ) and (fT7D, the SMC study 



of the semi-inclusive processes hints that the antiquarks in the sea are not 
strongly polarized. 

2.1.6 Baryon magnetic moments 

One of the puzzling aspects of the proton spin problem is that, given the 



significant deviation of the quark spin factors Ag in Eq.( 118 ) from the sQM 
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values, it is hard to see how could the same (Aq) sQM values manage to yield 
such a good description of the baryon magnetic moments, as shown in Table 
1. 

For this we can only give a partially satisfactory answer : If we assume 
that the anitquarks in the proton sea is not polarized A^ = 0, for which 
the SMC result Eq. ( |127|) gives some evidence (and it is also a prediction 
of the chiral quark model to be discussed in Sec. 3), we can directly use 
the Ag of Eq. (|118|) to evaluate the polarization difference: Aq = A q = Aq 
in Eq . (f49l) . We can then attempt a fit of the baryon magnetic moments in 
exactly the same way we had fit them by using ( Aq) as in Table 1. The 

V / sQM 

resultant fit, surprisingly, is equally good — in fact better, in the sense of 
lower x 2 fUl UMI • Namely, both the sQM Ag and experimental values of Aq 



can, rather miraculously, fit the same magnetic moment data. In this sense, 
the new spin structure poses no intrinsic contradiction with respect to the 
magnetic moment phenomenology. 

That it is possible to fit the same baryon magnetic moments with (Ag) ^ 

and (Ag) ^ is due to the fact that the baryon moment, such as Eq. (|50|) , is 
a sum of products fi B = J2 Ag /i g hence different (Ag) s can yield the same 

/ig if (/t g ) s are changed correspondingly. In both cases we have /i u = — 2fi d 
and fi s ~ F° r t ne sQM case, we find \i d ~ — 0.9/% while for the ex- 

perimental Ag case, we need /i d ~ — 1.4/%. This shift means a 35% change 
in the constituent quark mass value — thus a 35% difference with the con- 
stituent quark mass value obtained from the baryon mass fit in Eq.([39]). 
Consequently, we regard the magnetic moment problem still as an unsolved 
puzzle. 

2.2 DIS on proton vs neutron targets 
2.2.1 Lepton-nucleon scatterings 

The spin-averaged nucleon structure function F 2 can be expressed in terms 
of the quark densities as in Eq. (|108|) 



[x] = x 



F?(x) = x 



iJ {d + dj + liu + ^ + lis + s) 
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where we have used the isospin relations of (u) 
Their difference is 



(d)_ and (d) 



u 



where Z3 



x [Fi {X) 



(u — d) 



(u — d) + (u — d 
2X 3 + 2 (u - ctj 



u — d 



with it integral being the third component 

of the isospin: Jq dxX^ (x) = |. The simple assumption that u = d in the 
quark sea, which is consistent with it being created by the flavor-independent 
gluon emission, then leads the Gottfried sum rule [[29^ 



IG 



128) 



'OX o 

Experimentally, NMC found that, with a reasonable extrapolation in the 



very small-x region, the integral Ig deviated significantly from one-third |30 



T G = 0.235 ± 0.026 = \ + - , 

3 3 Jo 



u (x) — d (x) 



dx. 



(129) 



This translates into the statement that, in the proton quark sea, there are 
more d-quark pairs as compared to the u-quark pairs. 



u — d 



-0.147 ±0.026. 



(130) 



Remark Gottfried sum rule does not follow directly from QCD without ad- 
ditional assumption. Unlike the g\ sum rule, the Gottfried sum rule can 
not be derived from QCD via operator product expansion. A simple way 
to see this: Because the spin-independent structure function F 2 has op- 
posite crossing symmetry property from that of gi, only even-n terms 
can contribute. Hence there is no way to obtain a non-trivial relation 
for the odd-n moment sums of F2 (which the Gottfried sum rule would 
be an example). But in the context of parton model, the Gottfried sum 
provides us with an important measure of the flavor structure of the 
proton quark sea. 



2.2.2 Drell-Yan processes 

Because to conclude that NMC data showing a violation of the Gottfried 
sum rule one needs to make an extrapolation into the small-x regime, an 
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independent confirmation of u 7^ d would be helpful. A measurement of the 
difference of the Drell-Yan process of proton pN — > l + l~X on proton and 
neutron targets can detect the antiquark density because in such a process 
the massive pair is produced by (qq) annihilations |31 . 

Let us denote the differential cross sections as 

d 2 a {pN -> 1+rX) 



dy/rdy 

v q=u,a,s 



q p (xi) <? T (x 2 ) + <? P (»i) <? T (^2) 



(131) 



where a/t = ^= with -^/s being the CM collision energy and M is the invariant 
mass of the lepton pair, y being the rapidity, the fraction of momentum 
carried by the parton in the projectile (P) is given by x\ = \/re v and the 
fraction in the target (T) given by X2 = \fre~ v . Explicitly writing out the 
quark densities of Eq. ( |131| ) : 



a 



%ixa f 4 r . ,. . , . . 1 
— j- [u (xi) u (x 2 ) + u (xx) it (x 2 )] + - 



d (xi) d (x 2 ) + d (xi) d (x 2 ) 



a 



pi 1 



8ira J4 



9^ 19 



it (xi) d (x 2 ) + it (xi) 0? (x 2 ) 



1 r 
9 



d (xi) u (x 2 ) + ti (xi) it (x 2 ) 



+ s term 



+ s term 



In this way the DY cross section asymmetry can be found: 



A 



a 



pp 



a 



pn 



DY 



o-pp + cr pn 



[4m (xx) — d (xx)] 


u (x 2 ) - d (x 2 ) 


+ [it (x 2 ) - ci (x 2 )] 


4it (xi) — (i (xi) 


[4w (xi) + d(xi)] 


u (x 2 ) + d(x 2 ) 


+ [u (x 2 ) + d (x 2 )] 


4u (xi) + d (xi) 



(4A - 1) (A - l) + (A - 1) (4A - l) 
(4A + 1) (A + l) + (A + 1) (4A + l) 

where A (x) = it (x) jd (x) and A (x) 



(132) 



it (x) / d (x) . Thus with measurements 
of Ady and data fit for A in the range of (2.0, 2.7) , the NA51 Collaboration [[32 
obtained, at kinematic point of y 
antiquark distributions to be 



and Xi = x 2 



x 



0.18, the ratio of 



u/d = 0.51 ±0.04 ±0.05 (133) 

confirming that there are more (by a factor of 2) <i-quark pairs than it-quark 
pairs. 
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3 The Proton Spin-Flavor Structure in the 
Chiral Quark Model 

3.1 The naive quark sea 

A significant part of the nucleon structure study involves non-perturbative 
QCD. As the structure problem may be very complicated when viewed di- 
rectly in terms of the fundamental degrees of freedom (current quarks and 
gluons), it may well be useful to separate the problem into two stages. One 
first identifies the relevant degrees of freedom (DOF) in terms of which the 
description for such non-perturbative physics will be simple, intuitive and 
phenomenologically correct; at the next stage, one then elucidates the rela- 
tions between these non-perturbative DOFs in terms of the QCD quarks and 
gluons. Long before the advent of the modern gauge theory of strong inter- 
action, we have already gained insight into the nucleon structure with the 
simple nonrelativistic constituent quark model (sQM). This model pictures a 
nucleon as being a compound of three almost free u- and d- constituent quarks 
(with masses, much larger than those of current quarks, around a third of the 
nucleon mass) enclosed within some simple confining potential. There are 
many supporting evidence for this picture. We have reviewed some of this in 
Sec. 1. Also, the nucleon structure functions in the large momentum fraction 
x region, where the valence quarks are expected to be the dominant physical 
entities, are invariably found to be compatible with them being evolved from 
a low Q 2 regime described by sQM. For this aspect of the quark model we 



refer the reader to Ref. 33 



However in a number of instances where small x region can contribute 
one finds the observed phenomena to be significantly different from these 
sQM expectations. This has led many people to call sQM the "naive quark 
model" and to suggest a rethinking of the nucleon structure. But we would 
argue that the approach is correct, and only the generally expected features 
of the quark-sea are too simple. This "naive quark-sea" (nQS) is supposed 
to be composed exclusively of the u and d quark pairs. Namely, based on the 
notion of OZI rule, one would anticipate a negligibly small presence of the 
strange quark pairs inside the nucleon. This implies, as given in Eq.(p^), a 
pion-nucleon sigma term value of a^N — 25 MeV. Furthermore, the similarity 
of the u and d quark masses and the flavor-independent nature of the gluon 
couplings led some people to expect that d = u, thus to the validity of the 
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Gottfried sum rule, Eq.( 

In the sQM, there is no quark-sea and the proton spin is build up entirely 
by the valence quark spins. We have deduced the quark contributions to the 
proton spin as in Eq.([|7]), which leads to an axial- vector coupling strength 
of gA = Au — Ad = 5/3. If one introduces a quark-sea, the nQS feature of 
s ~ (thus As ~ 0) leads us to the Ellis- Jaffe sum rule, Jq 1 dxg\ (x) = 0.175. 

Phenomenologically none of these nQS features 

Features of the naive quark sea 
flavor : s = and d = u 

spin : As = (A ? = Ag) sea 

have been found to be in agreement with experimental observations. As far 
back as 1976, the connection of the ct^n value to the strange quark content 
of the nucleon has been noted. It was pointed out that the then generally 
accepted phenomenological value of 60MeV differed widely from the OZI 
expectation flOfl. In recent years, the a^N value has finally settled down to 
a more moderate value of a^N — 45 MeV when a more reliable calculation 



confirmed the existence of a significant correction due to the two-pion cut JTT] 



Nevertheless, this reduced value still translates into a nucleon strange quark 
fraction of 0.18, see Eq.(|7|). 

As for the proton spin, starting with EMC in the 1980's, the polarized 
DIS experiments of leptons on proton target have shown that Ellis- Jaffe sum 
rule is violated. The first moment the spin-dependent structure function gi 
has allowed us to obtain the individual Aq of Eq.( |118j ). We have already 
noted that they are all less than the sQM values of Eq . (f47l) > suggesting that 
for each flavor the quark-sea is polarized strongly in the opposite direction 
to the proton spin. 

A<? = (Aq) sQM + (Aq) sea <(Aq) sQM 



Furthermore, the recent SMC data on the semi-inclusive DIS scattering |26 
tentatively suggested A„ ~ Aj ~ 0. Thus while the inclusive experiments 
point to a negatively polarized quark sea, the semi-inclusive result indicates 
that the antiquarks in this sea are not polarized. 

The NMC measurement of the muon scatterings off proton and neutron 
targets shows that the Gottfried sum rule is violated JHJ. It has been inter- 
preted as showing d > u in the proton. This conclusion has been confirmed 
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by the asymmetry measurement (by NA51[32|) in the Drell-Yan processes 



with proton and neutron targets, which yield, at a specific quark momentum 
fraction value (x = 0.18), the result of d ~ 2u in Eq.( |133| ). 

To summarize, the quark- sea is "observed" to be very different from nQS. 
It has the following flavor and spin structures: 

Observed features of the quark sea 
flavor : d > u and s 7^ 

spin : (Aq) sea < yet Ag ~ 0. 

By the statement of s ^ 0, we mean that OZI rule is not operative for the 
strange quark. Recall our discussion in Sec. 1, this means that the couplings 
for the (ss)-pair production or annihilation are not suppressed, although the 
process may well be inhibited by phase space factors. Namely, a violation of 
the OZI rule implies that, to the extent one can ignore the effects of SU(3) 
breaking, there should be significant amount of (ss)-pairs in the proton. 

3.2 The chiral quark idea of Georgi and Manohar 

Let us start with theoretical attempts to understand the flavor asymmetry 
of d > u in the proton's quark sea: 

Pauli exclusion principle and the u-d valence-quark asymmetry in the 
proton would bring about a suppression of the gluonic production of v! s 
(versus d! s). Thus it has been pointed out long ago[B3J that d = u would not 



strictly hold even in perturbative QCD due to the fact the u's and d's in the 
qq pairs must be antisymmetrized with the u's and d's of the valence quarks. 
This mechanism is difficult to implement as the parton picture is intrinsically 
incoherent. In short, the observed large flavor- asymmetry reminds us once 
more that the study of quark sea is intrinsically a non-perturbative problem. 
Pion cloud mechanism ]35| is another idea to account for the observed 



d > u asymmetry. The suggestion is that the lepton probe also scatters off 



the pion cloud surrounding the target proton (the Sullivan process [[36|), and 
the quark composition of the pion cloud is thought to have more ds than u s. 
There is an excess of 7r + (hence d's) compared to 7r~, because p — > n+ii + , but 
not a 7T _ if the final states are restricted nucleons. (Of course, tt°s has d — u.) 
However, it is difficult to see why the long distance feature of the pion cloud 
surrounding the proton should have such a pronounced effect on the DIS 
processes, which should probe the interior of the proton, and also this effect 
should be significantly reduced by the emissions such as p — > A ++ + n~ , etc. 
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Nevertheless, we see that the pion cloud idea does offer the possibility to 
getting a significant d > u asymmetry. One can improve upon this approach 
by adopting the chiral quark idea of Georgi and Manohar pTfl so that there 
is such a mechanism operating in the interior of the hadron. Here a set 
of internal Goldstone bosons couple directly to the constituent quarks inside 
the proton. In the following, we will first review the chiral quark model which 
was invented to account for the successes of simple constituent quark model. 



The chiral quark idea Although we still cannot solve the non-perturbative 
QCD, we are confident it must have the features of (1) color confinement, 
and (2) spontaneous breaking of chiral symmetry. 

Confinement: Asymptotic freedom a s (Q) — > suggests that the run- 

ning coupling increases at low momentum-transfer and long distance, and 
a s (Aqcd) — 1 is responsible for the binding of quarks and gluons into 
hadrons. Experimental data indicates a confinement scale at 

Aqcd ^ 100 to 300 MeV. (134) 

Chiral symmetry breaking: There are three light quark flavors, m Uy d, s < 
h-QCD- In the approximation of m u ^, s — 0, the QCD Lagrangian is invari- 
ant under the independent SU(3) transformations of the left-handed and 
right-handed light-quark fields. Namely, the QCD Lagrangian has a global 
symmetry of SU(3)l x SU(3)r. If it is realized in the normal Wigner mode, 
we should expect a chirally degenerate particle spectrum: an octet of scalar 
mesons having approximately the same masses as the octet pseudoscalar 
mesons, spin ~~ baryon octet degenerate with the familiar | + baryon octet, 
etc. The absence of such degeneracy suggests that the symmetry must be 
realized in the Nambu-Goldstone mode: the QCD vacuum is not a chiral 
singlet and it possesses a set of quark condensate (0 \qq\ 0) ^ 0. Thus the 
symmetry is spontaneously broken 

SU{3) L x SU{3) R - SU(3) L+R 

giving rise to an octet of approximately massless pseudoscalar mesons, which 
have successfully been identified with the observed (tt, K, rj) mesons. 

The QCD Lagrangian is also invariant under the axial U(l) symmetry, 
which would imply the ninth GB m v i ~ m v . But the existence of axial 
anomaly breaks the symmetry and in this way the eta prime picks up an 
extra mass. 
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Both confinement and chiral symmetry breaking are non-perturbative 
QCD effects. However, they have different physical origin; hence, it's likely 
they have different distance scales. It is quite conceivable that as energy 
Q decreases, but before reaching the confinement scale, a s (Q) has already 
increased to a sufficient size that it triggers chiral symmetry breaking (xSB). 
This scenario 

A QCD <A xSB ~lGeV. (135) 

is what Georgi and Manohar have suggested to take place. The numerical 
value is a guesstimate from the applications of chiral perturbation theory: 
A x sb — 47r f v with f n being the pion decay constant. Because of this separa- 
tion of the two scales, in the interior of hadron, 

Aqod < Q < A-xSB, 

the Goldstone boson (GB) excitations already become relevant (we call them 
internal GBs), and the important effective DOFs are quarks, gluons and in- 
ternal GBs. In this energy range the quarks and GBs propagate in the QCD 
vacuum which is filled with the qq condensate: the interaction of a quark 
with the condensate will cause it to gain an extra mass of ~ 350 MeV. This 
is the chiral quark model explanation of the large constituent quark mass, 
(much in the same manner how all leptons and quark gain their Lagrangian 
masses in the standard electroweak theory). The precise relation between the 
internal and the physical GBs is yet to be understood. The non-perturbative 
strong gluonic color interactions are presumably responsible for all these ef- 
fects. But once the physical description is organized in terms of the resultant 
constituent quarks and internal GBs (in some sense, the most singular parts 
of the original gluonic color interaction) it is possible that the remanent inter- 
actions between the gluons and quarks/ GBs are not important. (The analogy 
is with quasiparticles in singular potential problems in ordinary quantum me- 
chanics.) Thus in our xQM description we shall ignore the gluonic degrees 
of freedom completely. 

Remark One may object to this omission of the gluonic DOF on ground 
that the one gluon exchange is needed to account for the spin- 
dependent contributions to the hadronic mass as discussed in Sec. 1. 
However, in the xQM the constituent quarks interact through the ex- 
change of GBs. The axial couplings of the GB-quark couplings reduce 
to the same s ' Sj effective terms as the gluonic exchange couplings. 
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For a more thorough discussion of hadron spectroscopy in such a chiral 
quark description see recent work by Glozman and Riska|38 . 



3.3 Flavor-spin structure of the nucleon 

In the chiral quark model the most important effective interactions in the 
hadron interior for Q < 1 GeV are the couplings of internal GBs to con- 
stituent quarks. The phenomenological success of this model requires that 
such interactions being feeble enough that perturbative description is appli- 
cable. This is so, even though the underlying phenomena of spontaneous 
chiral symmetry breaking and confinement are, obviously, non-perturbative. 



3.3.1 Chiral quark model with an octet of Goldstone bosons 



Bjorken[|39|, Eichten, Hinchliffe and Quigg^Oj are the first ones to point out 
that the observed flavor and spin structures of nucleon are suggestive of the 
chiral quark features. In this model the dominant process is the fluctuation 
of a valence quark q into quark q' plus a Goldstone boson, which in turn is 
a (qq') system: 

Q± — > GB + q' T — > (qq') Q q' T - (136) 

This basic interaction causes a modification of the spin content because a 
quark changes its helicity (as indicated by the subscripts) by emitting a 
spin-zero meson in P-wave. It causes a modification of the flavor content 
because the GB fluctuation, unlike gluon emission, is flavor dependent. 

In the absence of interactions, the proton is made up of two u quarks and 
one d quark. We now calculate the proton's flavor content after any one of 
these quarks turns into part of the quark sea by "disintegrating", via GB 
emissions, into a quark plus a quark- ant iquark pair. 

Suppressing all the space-time structure and only displaying the flavor 
content, the basic GB-quark interaction vertices are given by 



d = g 8 q$q = g 8 ( u d 



9h 



dix + sK + u 
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u + 
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1 u 


K° 




) 





(137) 



41 



Thus after one emission of the u quark wavefunction has the components 



dn + + sK + + u I + -4= 

\V2 Vg, 



(138) 



which can be expressed entirely in terms of quark contents by using 7r + = ud, 
and K + = us, etc. Since tt° and r\ have the same quark contents, we can add 
their amplitudes coherently so that 




2 

-uu 
3 



1„- 1 
-dd H — ss 
3 3 



(139) 



Square the wavefunction we the obtain the probability of the transitions: for 
example, 



Prob 



u + 



7T+d_ 



(lid) d. 



(140) 



which will be used to set the scale for other emissions. At this stage we shall 
assume SU(3) symmetry. Hence all processes have the same phase space, 
and are proportional to the same probability a oc |gs| 2 - The specific values 
are listed in the 3rd column of Table 2. The 2nd column is the isospin 
counter-part obtained by the exchange of u <-> d : 



u + — > 


d + — ► 


SU(3) sym prob 
octet GB 


broken U(3) prob 
nonet GB 


u + — > (wd) Q gL 
u + — > (ms) s_ 

M + — > (uu) Q 

u + — > (dd^J^U- 

U + — > (ss) M_ 


d + — > (du) m_ 
d + — * (ds) s_ 
d + — > (dd) d_ 
d + — > (mm) d_ 
d + — > (ss) d_ 


a 
a 

4 

9 a 

\a 


a 

e 2 a 

( £± f ±2 ) 2 « 
(^)\ 

(¥) 2 « 



Table 2 xQM transition probabilities calculated in models with an octet 
GB in the SU(3) symmetric limit and with nonet GB and broken-U(3) 

breakings. 



Flavor content calculation From Table 2, one can immediately read off 
the antiquark number q in the proton after one emission of GB by the initial 
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valence quarks (2u + d) in the proton: 



4 1 

u = 2 x -a + a + -a = 2a (141) 
9 9 K ' 

, 1 \ 4 8 
a = 2 x a H — a M — a = -a 

9/9 3 

1\ / 1 \ 10 
2x a H — a + a H — a = — a 



9 J V 9 / 3 

Since the quark and antiquark numbers must equal in the quark sea, we have 
the quark numbers in the proton: 

u = 2 + u, d=l + d, 8 = S. (142) 

Spin content calculation GB emission will flip the helicity of the quark 
as indicated in the basic process of ( |136| ), while the quark- antiquark pair 
produced through the GB channel are unpolarized: 

^ (GB) = A=ty (?+) i> (t) ~ i> (q-) i> (<?+)] ■ (143) 

One of the first xQM predictions about the spin structure is that, to the 
leading order, the antiquarks are not polarized: 

A,- = = 0. (144) 



Before GB emissions as in Ql36| ), the proton wavefunction is given by Eq. 
giving the spin-dependent quark numbers in Eq . (f45|) . Now from the 3rd 
column in Table 2, we can read off the first-order probabilities: 

2 

Pi (u+ — > gL) = a Pi (u+ — > s_) = a Pi (u+ — > u_) = -a, (145) 

3 

or write this in a more compact notation as 

Pi (u+ ^) = (d_ + s_ + ^_)a. (146) 

From this we can also immediately obtain the related probabilities of Pi (u_ — > 
Pi (d + — >) , and Pi (gL — >•) . The sum of the three terms in Eq. (|145 ) being |a, 



the probability of no GB emission must then be (1 — |a). Combining the 
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Oth and 1st order terms of Eqs.(f4"oD and ( |146| ), we find the spin-dependent 
quark densities (coefficients in front of q±): 

1 8 \ 1 1, 2 , \ 5,, 2 . 

3/ \3 U+ + 3 U ~ + 3 + + 3 ) + IT +s - + 3 M -) a 

1 2 1 2 2 2 

+ + s + + -u+)a + -(w_ + s_ + -cL)a + + s + + -d+)a 



Together with Eq.( 144 ), we can then calculate the quark polarization in the 



proton Aq = A q + A q = A q = q + — q_ : 

4 37 1 2 

Au = a, Ad= a, As = -a. (147) 

3 9 3 9 v ; 

In order to account for the NMC data of Eq. ( |130| ) by u — d = — |a as in 
Eq.( |141| ), we need a probability of a ~ 0.22. But such a large probability 
would lead to spin content description that can at best be described as fair. 
For example it give a negative- valued total quark value of AE = 1 — 16a/3 ~ 
—0.17, which is clearly incompatible with the current phenomenological val- 
ues in Eq. ( [L18|) — although it was still marginally consistent with the original 
EMC value when this calculation was first performed [40]. Also, the antiquark 
numbers in Eq.( |141| ) leads to a fixed ratio of u/d — 0.75, which is to be com- 
pared to the NA51 result of 0.51, as given in Eq. (|133|) . 



3.3.2 Chiral quark model with a nonet of Goldstone bosons 



We have proposed|4]J a broken-U(3) version of the chiral quark model with 
the inclusion of the ninth GB, the r/meson. 

Besides the phenomenological considerations discussed above, we have 
also been motivated to modify the original xQM by the following theoretical 
considerations. It is well-known that l/N co ; or expansion can provide us with 
a useful guide to study non-perturbative QCD. In the leading l/N coZor expan- 
sion (the planar diagrams), there are nine GBs with an U(3) symmetry. Thus 
from this view point we should include the ninth GB, the t]' meson. However 
we also know that if we stop at this order, some essential physics would have 
been missed: At the planar diagram level there is no axial anomaly and 7/ 
would have been a bona fide GB. Also, it has been noted by Eichten et al.|4"0[ 



that an unbroken U(3) symmetry would also lead to the phenomenologically 
unsatisfactory feature of a flavor-symmetric sea: u = d = s, which clearly vi- 
olates the experimental results of Eqs.(|130|) and (|133j). Mathematically, this 
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flavor independence comes about as follows. Equating the coupling constants 
g$ = gi in the vertex which generalizes the coupling in Eq. (|137|) 



8=1 



£i = 9s E Q^i&iV + \h9iQV'<l ( 148 ) 



(Xi4>i = $ with Aj being the Gell-Mann matrices) and squaring the amplitude, 
one obtains the probability distribution of 

8 2 

E + o (ft) (??) ( 149 ) 

t=i 6 

which has the index structure as 

8 2 

E Mat + = (150) 

where we have use a well-known identity of the Gell-Mann matrices to obtain 
the equality. This clearly shows the flavor independence nature of the result. 

Calculation in the degenerate mass limit All this shows that we should 
include the ninth GB but, at the same time, it is crucial that this resultant 



flavor-U(3) symmetry be broken. In our earlier publication [41] we have im- 
plemented this breaking in the simplest possible manner by simply allowing 
the octet and singlet couplings be different. Namely, in the first round cal- 
culation, we stayed with approximation of m n = m s and a degenerate octet 
GBs. In this way we were able to show that with a choice of 

C= — ^ -1 (151) 

#8 

this broken U(3) xQM can account for much of the observed spin and flavor 
structure, see Column-5 in Table 3. 

Our calculation has been performed in the SU(3) symmetric limit (i.e. 
assumed all phase space factors are the same). In this spirit we have chosen 
to work with \gi\ = \g&\ . The relative negative sign is required primarily to 
yield an antiquark relation of d ~ 2u : as the model calculation gives a ratio 

u/d= C ' + * + 6 . (152) 
1 C +8 K J 
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Therefore, the experimental value of Eq.( |133| ) implies a negative coupling 
ratio : —4.3 < C < —0.7. We remark that the relative sign of the couplings is 
physically relevant because of the interference effects when we coherently add 
the rj contribution to those by rj and 7r°. After fixing this ratio, there is only 
one parameter a that we can adjust to yield a good fit. It is gratifying that 
a = 0.11 is indeed small, fulfilling our hope that once the singular features 
of the nonperturbative phenomenon of spontaneous symmetry breaking are 
collected in the GB degrees of freedom, the remanent dynamics among these 
particles is perturbative in nature. 

It should also be noted that we have compared these SU(3) symmetric 
results to phenomenological values which have been extracted after using 
the SU(3) symmetry relations as well. For example the result in Eq.( |118| ) 
have been extracted after using the SU(3) symmetric F/D ratio for hyperon 
decays as in Eq.( |115[ ). Similarly, we obtained a strange quark fraction value 
F (s) ~ 0.19 very close to that given in Eq. (|79| ) which was deduced from a n N 
and an SU(3) symmetric F/D ratio for baryon masses, Eq.([77]). Agreements 
are in the 20% to 30% range, indicating that the broken-U(3) chiral picture 
is, perhaps, on the right track. 



SU(3) and axial-U(l) breaking effects The quark mass difference m s > 
m n , and thus the GB nondegeneracy , would affect the phase space fac- 
tors for various GB emission processes. Such SU(3) breaking effects will 



be introduced [^5| (|6| in the amplitudes for GB emissions, simply through 
the insertion of suppression factors: e for kaons, 5 for eta, and ( for eta 
prime mesons, as these strange quark bearing GB's are more massive than 
the pions. Thus the probability a oc |g§| 2 are modifies for processes involving 
strange quarks, as shown in the last column of Table 2. The suppression 
factors enter into the probabilities for u + — > (uu) Q u_ and u + — > (dd^j^U- 
processes, etc. because they also receive contributions from the 77 and rj 



GBs. Following the same steps as those in Eqs. (|137| ) to ( |140| ), we obtain the 
probabilities as listed in the 4th column of Table 2. In this way the following 
results are calculated: 



d 



12 
1 

12 

1 r 



(2C + 5 + l) 2 + 20 
(2C + 5- l) 2 + 32 
(C - S) 2 + 9e 2 l a. 



(153) 
(154) 
(155) 
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and 



Au 

Ad 
As 



4 21 + A5 2 + 8C 2 + 12e 2 
3 9 

6 - S 2 - 2C 2 - 3e 2 



1 

~3 ~ 
— e 2 a 



(156) 

(157) 
(158) 



In the limit of ( = (i.e. no rj') and e = 5 = 1 (no suppression in the 
degenerate mass limit) these results are reduced to those of Eqs.( |141| ) and 

Results of the numerical calculation are given in the last column in Table 
3. Again our purpose is not so much as finding the precise best-fit values, but 
using some simple choice of parameters to illustrate the structure of chiral 



quark model. For more detail of the parameter choice, see Ref.[46 . 



Phenomenological 
value 



Eq. 

# 



Naive 
QM 



a = 



XQM 
SU 3 sym 
e = 6 = 

-C = i 

a = 0.11 



brok'n SU3 

t = 5 = 
-2C = 0.6 
a = 0.15 



u — d 


0.147 ±0.026 


u/d 


(0.51 ±0.09) x=0 . 18 


2s/ (u + d) 


~ 0.5 


cr nN : F (s) 


0.18 ±0.06(|?) 


F(3) /F(8) 


0.23 ±0.05 


9a 


1.257 ±0.03 


1 / / axial 


0.575 ±0.016 


(3F - D)a 


0.60 ±0.07 (|?) 


Au 


0.82 ±0.06 


Ad 


-0.44 ±0.06 


As 


-0.11 ±0.06 (|?) 


Au, Ad 


-0.02 ± (.11) 



a) 

1331) 



0? 
1? 
0? 



0? 
1 

3 



4 

3 



0.146 
0.56 
1.86 



0.19 

1 
3 

1.12 

2 
3 

0.67 



0.78 
-0.33 
-0.11 







0.15 
0.63 
0.60 

0.09 
0.22 

1.25 
0.57 
0.59 

0.85 
-0.40 
-0.07 
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Table 3 Comparison of xQM with phenomenological values. The 3rd 
column gives the Eq. numbers where these values are discussed. From there 
one can also look up the reference for the source of these values. Possible 
downward revision of the results by SU(3) breaking effects, as discussed in 
the text, are indicated by the symbol (J,?) . Those values with a question 
mark (?) in the 4th column are not strictly the sQM predictions, but are the 
common expectations of, what has been termed in Sec. 3.1, the "naive quark 



sea" . 



Since a SU(3) symmetric calculation would not alter the relative strength 
of quantities belonging to the same SU(3) multiplet, our symmetric cal- 
culation cannot be expected to improve on the naive quark model, i.e. 
SU(6), results such as the axial vector coupling ratio F/D = 2/3, which 
differs significantly from the generally quoted phenomenological value of 
F/D = 0.575 ± 0.016. To account for this difference we must include the 
SU(3) breaking terms: 

F Au - As 



D Au + As - 2 Ad 

2 6-a(25 2 + 4C 2 + |( 3e2 + 21 : 

3 6 -a(25 2 + 4C 2 + 9e 2 + 3) 



(159) 



Similarly discussion holds for the F j D ratio for the octet baryon masses. Here 
we choose to express this in terms of the quark flavor fractions as defined by 
Eqs.PD and (M): 



F (3) F (u) — F (d) 1 + 2 (u - d 



F(8) F(u) + F(d)-2F(s) 3 + 2 (u + d - 2s 
1 3 + 2a [2( + 5 - 3] 



3 3 + 2a[2C5+|(9-5 2 -12e 2 



(160) 



In the SU(3) symmetry limit of 5 = e = 1, we can easily check that Eqs.( |159| ) 
and (|160|) reduce to their naive quark model i.e. SU(6) values, independent 
of a and (. Again it is gratifying to SGG, clS displayed in Table 3, that xQM 
has just the right structure so the SU(3) breaking modifications make the 
correction in the right direction. 
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3.4 Strange quark content of the nucleon 

We have already discussed the number s of strange quarks in the nucleon 
quark sea and their polarization As. They are examples of the proton matrix 
elements of operators bilinear in the strange quark fields (p\sTis\p) , or in 
general we need to study the quark bilinear matrix elements of (p \qTiq\p) : 



3.4.1 The scalar channel 



This operator counts the number of quarks plus the number of antiquarks 



in the proton. In particular the octet components of (p 



uu 



dd 



and 



V 



uu + dd — 2ss 



have shown in Eq. 



p ) can be gotten by SU(3) baryon mass relations as we 
. But in order to separate out the individual terms, say 



uu + dd + ss 



P 



This is 



{p \ ss\ p), we would need the singlet combination [p 
provided by a n N which is a linear combination of the singlet and octet pieces. 
That is why a measurement of u^n allows us to do an SU(3) symmetric 
calculation of the strange quark content of the nucleon. 

We have emphasized that OZI violation means that the couplings for ss 
pair creation and annihilation may not be suppressed even though the phase 
space surely does not favor such processes. But the phase space suppression 
is a "trivial" SU(3) breaking effect. Our chiral quark model calculation is a 
concrete realization of this possibility: Had we ignored the phase space dif- 
ference, the GB-quark couplings are such that there would be more strange 
quark pairs than either of the nonstrange pairs in the quark sea, as ss pro- 
duction by either u or d valence quarks are not disfavored. Thus Eq. (|141|) 
give a relative quark abundance in the quark sea of 



u : d : s = 3 : 4 : 5 



(161) 



In the physical quark sea we do not really expect strange quark pairs to 
dominate because of their production is suppressed by SU(3) breaking effects. 

The xQM naturally suggests that the nucleon strange quark content s and 
polarization As magnitude are lowered by the SU(3) breaking effects as they 
are directly proportional to the amplitude suppression factors, see Eqs. ( |155| ) 
and ( |158| ). This is just the trend found in the extracted phenomenological 
values. GasserETI, for instance, using a chiral loop model to calculate the 



SU(3) breaking correction to the Gell-Mann-Okubo baryon mass formula, 
finds that the no-strange-quark limit-value of (ova^o * s modified from 25 to 
35MeV, [i.e. the baryon mass M 8 in Eq . (ffTD changed from —200 by SU(3) 
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breakings to —280 MeV], thus the fraction F (s) from 0.18 to 0.10. It matches 
closely our numerical calculation with the illustrative parameters, see Table 
3. 

The strange quark content can also be expressed as the relative abundance 
of the strange to non-strange quarks in the sea, which in this model is given 

as 

X s = . - x = 4 (C ~^ )2 2 +9£2 ~ 1.6e 2 = 0.6. (162) 
l(u + d) (2C + 5) 2 + 27 

This can be compared to the strange quark content as measured by the CCFR 
Collaboration in their neutrino charm production experiment |}4"8"| 



(xs) 

h 



0.477±0.063, where (xq) = / xq(x)dx, (163) 

Jo 



which is often used in the global QCD reconstruction of parton distributions [49 



The same experiment found no significant difference in the shapes of the 



strange and non-strange quark distributions f48 

[xs (x)] oc (1 — x) a 



xu (x) + xd (x) 



with the shape parameter being consistent with zero, a = — 0.02±0.08. Thus, 
it is reasonable to use the CCFR findings to yield 

A s ~ « ~ (164) 

which is a bit less than, but still compatible with, the value in Eq. (|162[) . 

Thus it is seen that the xQM can yield a consistent account of the strange 
quark content s of the proton sea. SU(3) breaking is the key in reconciling the 
s value as measured in the neutrino charm production and that as deduced 
from the pion nucleon sigma term. 

3.4.2 The axial-vector channel 

This operator measures the quark contribution to the proton spin. In par- 
ticular the octet components of 



p,s 



p,s 
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and 

(p, s 1^7^75" + ^7 M 7 5 d - 2 ^75 s | P, s 

can be gotten by SU(3) relations among the axial vector couplings of octet 
baryon weak decays, Eqs.( |114Q and ( |115[ ). But in order to separate out the 



individual terms, say (p, s s / y fl / y 5 s p, sy = 2s^As, we would need the singlet 
combination Au + Ad + As. This is provided by the first-moment of the 
structure function / g\dx which is a linear combination of the singlet and 
octet pieces. That is why a measurement of g\ (x) allows us to do an SU(3) 
symmetric calculation of the strange quark content of the nucleon. 

A number of authors have pointed out that phenomenologically extracted 
value of strange quark polarization As is sensitive to possible SU(3) breaking 
corrections. While the effect is model-dependent, various investigations [[42] - 



|H all conclude that SU(3) breaking correction tends to lower the magnitude 
of As. Some even suggested the possibility of As ~ being consistent with 
experimental data. Our calculation indicates that, while As may be smaller 
than 0.10, it is not likely to be significantly smaller than 0.05. To verify this 
prediction, it is then important to pursue other phenomenological methods 
that allow the extraction of As without the need of SU(3) relations. 

Besides polarized DIS of charged lepton off nucleon, we can also use other 
processes to determine As. In elastic neutrino-proton scattering, we can sep- 
arate out the axial form factors at zero momentum transfer, 



P' 



97u7 5 9 



P 



2u (p') 



G[ q) (Q 2 ) 7,75 + (Q 2 



u (p) . (165) 



Thus we have Gi (0) = Aq. The axial vector matrix element arises from 
Z-boson exchange is proportion to 



P 



p 



i 



p 



P 



(166) 



S7 M 7s s 



P 



where T 3 is the 3rd component of the weak isospin operator. The <^p 
can be separated out because the first two terms are fixed by the neutron 
axial coupling g&. Present data still have large error, however they are con- 
sistent with a As ^ 0J50]|. 

The measurements of longitudinal polarization of A in the semi-inclusive 
process of vN — > fiA + X [[H]] have also given support to a nonvanishing and 
negative As. In this connection, it's also important to pursue experimen- 
tal measurements to check the xQM prediction for a vanishing longitudinal 
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polarization of A in the semi-inclusive processes reflecting the proton spin 
property of As = 0. 



3.4.3 The pseudoscalar channel 

The nucleon matrix elements of the pseudoscalar quark density may be phys- 
ically relevant in Higgs coupling to the nucleon|^2[, etc. Such operators may 



be related to the axial vector current operator through the (anomalous) di- 
vergence equation ( |121| )p5|. If we define 



P 



(p\qh 5 q\p) = v q u{p)ij 5 u(p) (167) 
p) = -Ag2M p u(p)i l5 u(p) (168) 



trG"^ 



so that the non-strange divergence equations may be written as 



a. 



2M p Au = 2m u v u - 2M p [^ A 9J ( 16 9) 
2M p Ad = 2m d p d -2M p (^Ag\ 

We would need one more condition in order to separate out the individual 
m q v q terms. This may be obtained by saturation of the nonsinglet channel by 
Goldstone poles. Let us recall that the Goldberger-Treiman relation can be 
derived in the charge channel by the 7r =l= pole-dominance of the pseudoscalar 
density. After taking the nucleon matrix element of 



9^ {ui p l 5 dj = (m u + m d ) (m7 5 d) 

one obtains 

2M p g A = 2f w g wNN + (i ± (170) 

where fi± denotes the correction to the 7r =1= pole-dominance, and is the cor- 
rection to the gA as given by the GT relation. Repeating the same for the 
neutral isovector channel 

(ui p l 5 u - dj^sd) = 2m u (uij 5 u) - 2m d (dij 5 d) 

we have 

2M p g A = 2f 7T g 7rNN + H + {m u - m d ) {y u + v d ) (171) 

Comparing these two expressions for g& one concludes that the singlet den- 
sity iy u + u d ) must be small, on the order of correction to the GT expression 
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of qa- Assume that fi ~ // ± , thus v u = we can solve the two equations 
in ( |169|) in terms of the measured Am and Ad given in Eq.( |118| ). 



m u v u = 423 MeV m d v d = -761 MeV —Ag 

With these values we can also obtain 

m sVs = -451 MeV . 



-0.37. 



(172) 



(173) 



Because of the large strange quark masses m s , this translates into fairly small 
strange pseudoscalar matrix element of 



{p \si^ 5 s\ p) ~ 0.03 (p di^ h d 



P 



3.4.4 The vector channel 

Of course the vector charges 



Q' 



d?xV* (x) 



with VI 



~2 



Tly.—Q. 



(174) 



are simply the generators of the flavor SU(3). In terms of the form factors 
defined as 



P 



~2 



P 



u(p) 



(Q") 



u(p) 
(175) 



where Q 2 = {p' — p) 2 is the momentum transfer, we note that (^p 
are constrained by the quantum numbers of the proton: 



F} u) (0) - Ff* (0) = 1, Fl s, (0) = 0. 



<d) 



p 



(176) 



However, the magnetic moment form factor (0) needs not vanish. It 

is therefore interesting to measure this quantity. This can be done through 
the observation of parity violation in the scattering of charged-leptons off 
nucleon. The interference of the photon-exchange and Z-boson-exchange dia- 
grams can be used to isolate F^ (0) . For detailed discussion, see Refs. |54j] |5l| . 
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3.5 Discussion 



In these lectures we have described an attempt to understand the nucleon 
spin-flavor structure in the framework of a broken-U(3) chiral quark model. 
The broad agreement obtained with simple schematic calculations, as dis- 
played in Table 3, has been quite encouraging. If this approach turns out 
to be right, it just means that the familiar non-relativistic constituent quark 
model is basically correct — it only needs to be supplemented by a quark 
sea generated by the valence quarks through their internal GB emissions. 

Because the couplings between GB and constituent quarks are not strong, 
we can again use perturbation theory based on these non-perturbative de- 
grees of freedom — even though the phenomena we are describing are non- 
perturbative in terms of QCD Lagrangian quarks and gluons. Features such 
as d ~ 2u are seen to be clear examples of nonperturbative QCD physics, as 
they are quite inexplicable in terms of a quark sea generated by perturbative 
gluon emissions. (If one gets beyond the perturbative gluonic picture, this 
d 7^ u property is not peculiar at all, as the nucleon is not an isospin sin- 
glet and there is no reason to expect that its quark sea should be an isospin 
singlet.) 

In the case of the proton spin structure, because the most often discussed 
theoretical interpretation is the possibility of a hidden gluonic contribution, 
it has led some to think that other approaches, such as xQM, must be 
irrelevant. But the alternative theories are attempting a different description 
by using different degrees of freedom. To be sure, the QCD quarks and gluons 
are the most fundamental DOF. But we cannot insist on using them for such 
non-perturbative problems as the hadron structure. An analogy with the 
nucleon mass problem will illustrate our point. 

The canonical approach to study the various quarks/gluon contributions 
to the nucleon mass is through the energy-momentum trace anomaly equation 



9£ = m u uu + m d dd + m s ss - (11 - -n f ) ^trG^G^. (177) 

\ 6 / on ' 

Just like the more familiar axial vector anomaly equation, the naive diver- 
gence is given by quark masses while the anomaly term is given by the gluon 
field tensor. (Of course, here we are using the Lagrangian quark and gluon 
fields.) When taken between the proton states, this equations yields 

M p = m n (p uu + dd + m s (p\ss\p) + gluon term (178) 
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The first term is just the a^N — 45 MeV representing a tiny contribution by 
the nonstrange quarks, while the second term can also be estimated ||52||: 



rrisa^N 

m s {p \ss\p) 



2m. 



n 



3m n M 8 



250 MeV. (179) 



[Because this is an SU(3) calculation, the strange quark term is somewhat 
overestimated.] One way or other, we see that most of the proton mass came 
from the gluon term|p7||. Heavy quark terms can also be included but their 
contributions as explicit quark terms just cancel the corresponding heavy 
quark loops in the gluon terms. In this sense they decouple ||58||. 

This led to an important insight: nucleon mass is mostly gluonic. But 
in terms of the QCD quarks and gluons, it is difficult to say anything more. 
That is why the description provided by the constituent quark model is so 
important. In this picture much more details can be constructed: hyperfine 
splitting, magnetic moments, etc. 

The important point is that these two approaches are not mutually ex- 
clusive. While the constituent quark model does not refer explicitly to gluon, 
the above discussion suggests that it is the non-perturbative gluonic inter- 
action that brings about the large constituent quark masses. (In the xQM 
this takes the form of quark interaction with the chiral condensate of the 
QCD vacuum.) We believe that this complimentarity of the QCD and sQM 
descriptions holds for the flavor-spin structure problem as well. The non- 
perturbative features can be described much more succinctly if we use the 
non-perturbative DOF of constituent quarks and internal GBs. Thus it is 
quite possible that the statement of a significant gluonic contribution to the 
proton spin and a correct description of spin structure by the xQM can both 
be valid — just the same physics expressed in two different languages. 
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